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Abstract 

It is shown theoretically that after light storing in a medium of four-level 
atoms it is possible to release a new pulse of a different frequency, the process 
being steered by another driving beam. It is also possible to store one pulse 
and to release two different ones, with their time separation and heights being 
controlled. 
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It has been shown both theoretically and experimentally (see, e.g., Refs |L]-|5|] and a 
review ||) that a light pulse, propagating in a medium composed of three-level atoms in 
a A configuration, suitably driven by another pulse, can be stopped and later released 
in a controlled way. The process is interpreted in terms of inducing a transient Raman 
coherence between two lower atomic states or, in the language of quasi-particles, in terms 
of an adiabatic evolution of the so-called dark-state polariton . Natural questions arise 
whether it is possible to convert the stored light in a controlled way into a pulse of a frequency 
different from that of the stopped one, or into more pulses of different frequencies, or what 
happens if the Raman coherence is due to a transient absorption of more than one pulse. 
Some experimental results concerning light frequency conversion during the process of its 
storing and retrieving in a four-level system have recently been mentioned by Matsko et al. 
|J. In this paper we present a complete theoretical analysis of the situations when a four- 
level-atom medium is driven by two control laser fields and, after stopping one or two probe 
pulses, one or two pulses can be released. We also point out the possibility of controlling 
the process by a relative time shift of the control fields. 

Consider a quasi one-dimensional medium of four-level atoms with two lower metastable 
states | b > and |c > and two upper states \a > and \d > (Fig.0). The position of an atom 
is described by the variable z, which is considered continuous. The weak signal fields 1 and 
3 couple the initial state \b > with \a > and \d >, while the control fields 2 and 4 couple 
|c > with | a > and \d >, respectively. The interaction Hamiltonian is V — —dJ2^=i £j cos (ft j, 
with <pj = ujjt — kjZ, 6j = €j(z, t) being slowly varying envelopes and all the fields having the 
same linear polarization. The matrix elements of the dipole moment d\ = (d) ab , d 2 = {d) ac , 
d 3 = (d)bd, d>4 = (d)cd are taken real. Resonant conditions concerning all the couplings are 
assumed, i.e. huj\ = E a — E b , Tiu 2 = E a — E c , fuo 3 = Ed — E b , %uj 4 = Ed — E c . 

The evolution equation ihp = [H, p] for the density matrix p = p(z, t) for an atom 
at position z, after making the rotating-wave approximation, transforming-off the rapidly 
oscillating factors: p ab = a ab exp(-z^i), p ac = a ac exp(-z0 2 ), Pbc = cr bc exp[i(0i - (f) 2 )}, 
Pdb = cr db exp(-i(f) 3 ), p dc = cr (fc exp(-i0 4 ), p ad = a ad exp[i((f) 3 - fa)], p u = an, and after 
adding relaxation terms describing the spontaneous emission within the system, takes the 
form 

ih°aa = -~eidi(a ba - a ab ) - ^e 2 d 2 {cr ca - a ac ) - %%(Yl + T a c )a aa , 
ih& bb = --eidi(a ab - a ba ) - ^e 3 d 3 (a db - cr bd ) + ihTla aa + ihT^a d d, 
iha cc = -^e 2 d 2 {cr ac - a ca ) - ^d A (a d c - cr cd ) + ihT a c a aa + ihTfa dd , 

iho dd = -~e 3 d 3 (o- bd - a d b) - ^d 4 (a cd - a dc ) - iTi^t + ^t) a dd, 

1 1 1 ih 

ih& ab = --e x d x {a bb - a aa ) - -e 2 d 2 a cb + -e 3 d 3 a ad - — {VI + T a c )a ab , (1) 

ih&ac = --ei^i^fec - -^d 2 (a cc - a aa ) + -e 4 d 4 a ad - — (T^ + T a c )a ac , 

1 1 1 1 ih 

ifrcrad = -^\d\o bd - -e 2 d 2 a cd + -e 3 d 3 a ab + -t 4 d 4 a ac - y (Tt + T a c + T b + T c )a ad , 

1 . 1 j 1 j 1 , 

iha bc = --exdiCJac + -e 2 a 2 a ba - -e 3 d 3 a dc + -e 4 d 4 a bd , 
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incTbd = --£id\(Tad - -^d 3 (a d d - + -e4«4C& c - y (r 6 + r c )cr 6d , 

111 z?i 
^o"crf = —^■T.^od + -^d 3 a cb - -e^d 4 (a d d - o- cc ) - y (Tb + F c )<J c d, 

where T% is the decay rate of the state \a > to \b > etc. 

The propagation equations for the signal fields 1 and 3 are written as usual in the slowly 
varying envelope approximation || after rejecting the second space and time derivatives of 
€j. In the conditions of the resonance they read 

~R 1 ET = l N dl (T ab , 

oz c at tQC 
de 3 lde 3 u 3 

UZ C Ui tQC 

where N is the atom density, eo is the vacuum electric permittivity and use has been made of 
the fact that in the resonance conditions a ao and a d b are imaginary numbers. Similarly as in 
earlier papers, we have neglected propagation effects for the driving fields, i.e. e 2)4 = e 2j 4(i). 

Eqs (1) and (2) have been solved numerically in the moving window frame of reference: 
if = t — z j c, z' = z using the method described by Shore 0. Switching the driving fields on 
and/or off was modeled by a hyperbolic tangent. The initial probe pulse was taken as the 
sine square shape ei(0,i) = ei sin 2 [7r(t — t 2 )/(t2 — £j.)]©(t — ti)Q(t 2 — t), while the initial 
condition for the atomic part was <7bb{z, 0) = 1, with other matrix elements equal to zero. 

We have performed simulations for somewhat arbitrarily chosen data, being however of 
realistic orders of magnitude. The model atomic energies were E a = —0.10 a.u., Eb = —0.20 
a.u., E c = —0.18 a.u., E d = —0.05 a.u., all the relaxation rates due to the spontaneous 
emission were taken Y a b ' a c = 2.4 x 10" 9 a.u., from which the corresponding transition dipole 
moments were calculated. The length of the atomic sample was 3 x 10 7 a.u. (1.6 mm) and 
its density 3 x 10" 13 a.u. (2 x 10 12 cm 3 ). The initial signal pulse length was 10 11 a.u. (2.4 
/is) and eio = 10~ 10 a.u. (which corresponds to the power density of 3.5 xl0~ 4 Wcm" 2 ); the 
maximum value of the control field amplitudes was 1.2 x 10 -9 a.u. (50 mWcm -2 ). 

The upper part of Fig. |2| shows for comparison the results obtained by solving Eqs (1) 
and (2) in the case corresponding to the recently studied light storage in a single A system 
(cf. the results of Ref. M). The pulse's electric field is shown as a function of the "local" 
time t' (which for our data is almost equal to t). The left peak of the signal pulse is the 
untrapped fraction of the incoming pulse, i.e. its fraction transmitted by the medium before 
the control pulse has been switched off. The rest of the pulse is trapped and later released 
(the right peak) due to turning the control pulse back on. We have checked that, as expected, 
the time interval between the left and right peaks is the same as that between the instants 
of switching the control pulse off and on. 

The lower part of Fig. ^ shows the new result predicted in this paper. As before, the 
signal pulse e\ is trapped due to switching the control field e 2 off and its untrapped part is 
observed. However, in this case, after some time the other control field e 4 is turned on, which 
results in generating a new pulse e 3 of frequency different from that of the original trapped 
signal pulse. Also in this situation the instant of appearing of the released pulse (now the 
field 3) is controlled by choosing the moment of the switch-on of the driving field (now the 
field 4). Numerical calculations show that the height and shape of this new released pulse 
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strongly depend on the process parameters, in particular on the maximum value and slope 
of the control field 4. 

From our observations it follows that it should be possible to trap a single pulse and 
then to release two or more pulses. Such a possibility is demonstrated in Figs 3a and 
3b, which show the released pulses in the cases in which the driving field was switched 
on, respectively, before and after the driving field e-i. The difference of the pulses heights 
shows that by choosing the order of the turn-on of the driving pulses and their time delay 
we can influence the fraction of the atomic coherence taken over by each of the released 
pulses. We have also obtained the results (not shown here) which prove that it is possible 
to simultaneously store two pulses and release one or two new ones. 

Extrapolating the interpretation concerning a single A system one may look, analogously 
as in Ref. for the solutions of Eqs (1) and (2), assuming a perturbative, adiabatic and 
relaxationless evolution. It follows then that in particular 

a — _fl^l — e 3^3 ^ 

Those approximations allow one to find a shape-preserving solution of the corresponding 
Maxwell-Bloch equations ^(z,t) = — jQv(t')dt',t = 0), ^ being a combination of field 
and atomic variables 

^2^2 r 2NTlu)\ _i_ rf4£4Ui 

yj, _ di 1 eg bc d 3 UJ 3 3 

d\ 2 I 2Nhun , d^un 9 



The velocity v is given by 



d% c 2 1 1 d\^± 

d? 6 2 + +d^ 6 4 



V = C 1 3 (5) 

However, after switching the field 4 on, the solution \I/ tends to y^ e 3 instead of e 3 , which 
means that this solution cannot fully characterize the pulse release phase of the evolution. 

The latter would be described by an adiabatic evolution of ^' = This means that 

V ui 

light storing and releasing in a double A system must include a nonadiabatic phase; the 
measure of nonadiabaticity is the ratio of the frequencies of the two signals. This has been 
checked numerically: we have compared the computed coherence <Jb c (t') with its adiabatic 
approximations given by both parts of Eq. (3). While the evolution of o^c in the storage 
phase was well reproduced by the part of Eq. (3) including the fields 1 and 2, this was 
not the case for the part including the fields 3 and 4 in the release phase. The coherence 
calculated from the two parts of Eq. (3) exhibited a discontinuity, which could be reduced 
by artificially introducing the factor J^- The failure of the adiabatic approximation in 
this case is connected with the fact that a four level atom in the conditions of resonance, 
dressed by the four interactions (1-4) satisfying the second equality in Eq. (3), has a double 
real eigenvalue (equal to the bare energies). Thus the assumptions of the adiabatic theorem 
are not satisfied. Considerations analogous to those presented above but performed in the 
formalism of second quantization would allow an interpretation of the process in terms of 
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quasiparticles. As in the case of a single A both the storage and release phases can be seen 
as an adiabatic evolution of a dark state polaritons, but here the two polaritons in the two 
phases are not identical and a nonadiabatic transformation of one into another must take 
place. 

The above results suggest new interesting possibilities of controlling light propagation 
effects. A systematic quantitative analysis of the dependence of the pulse shapes, as well as 
of the time evolution of the atomic properties, on the numerous accessible control parameters 
will be the subject of a future work. 
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FIGURES 



FIG. 1. The double A scheme of levels and couplings; the indices 1 and 3 refer to signal fields 
and 2 and 4 - to control fields. 

FIG. 2. The shape of the signal pulses at the end of the sample as the function of the "local" 
time t'; also shown are the control pulses (for pictorial reasons their values were reduced by the 
factor of 0.03); (a) the case of a single A system: 63,4 = 0, e%, solid line; €2, dashed line; for 
comparison is also shown the shape of initial pulse, dotted line; (b) the case of a double A system: 
the transmitted part of the original signal pulse ei, solid line; the newly created signal pulse €3, 
dashed line; 62, short-dashed line; £4, dotted line. 

FIG. 3. The shape of the released pulses depending on the order of switching on of the control 
pulses: (a) the pulse (4) switched on before the pulse (2); (b) the pulse (4) switched on after the 
pulse (2). The line styles as in Fig. |2|b. 
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